range corresponding to turbulent inertial scales and explore how the turbulent forcing experienced 
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(N ; Abstract 

We use an acoustic Lagrangian tracking technique, particularly adapted to measurements in 
open flows, and a versatile material particles generator (in the form of soap bubbles with ad- 
justable size and density) to characterize Lagrangian statistics of finite sized, neutrally bouyant, 
particles transported in an isotropic turbulent flow of air. We vary the size of the particles in a 
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by the particles depends on their size. We show that, while the global shape of the intermittent 
O ' acceleration probability density function does not depend significantly on particle size, the accel- 
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£>y eration variance of the particles decreases as they become larger in agreement with the classical 

Qh. k 7 > 3 scaling for the spectrum of Eulerian pressure fluctuations in the carrier flow. 
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Particle laden flows play an important role in several industrial applications (diesel en- 
gines, mixers, separative techniques, etc.), environment issues (atmospheric dispersion of 
pollutants, sediments transport, etc.) and many natural systems (plankton in seas, spores 
and pollen dispersion, water droplets in clouds, planetesimals in accretion disks, etc.). In 
most cases the carrier phase is a turbulent flow with eddies distributed over a wide range 
of relevant scales (called the inertial range of turbulence) : from the energy injection scale 
(L) down to the Kolmogorov dissipation scale (rj). If the advected particles are neutrally 
buoyant and small (typically comparable in size or smaller than 77) they behave as tracers 
for fluid particles. Such tracers are commonly used for single phase hydrodynamic mea- 
surements such as PIV, LDV and PTV. However, in many practical and natural situations, 
the transported particles have a density different from that of the surrounding fluid and/or 
a size larger than 77. Their dynamics, which is then affected by particles inertia, deviates 
from that of fluid particles. The statistical description of inertial particles interacting with 
a turbulent flow remains a largely open question which clearly requires further experimen- 
tal insight. From a theoretical point of view, even writing (not to mention solving) the 
particle motion equation in the most general case remains a challenge and only limit cases 
have been treated at present. For instance, in most analytical approaches, particles are as- 
sumed to be small compared to the scales of turbulence, so that the fluid around is assumed 
locally laminar. In this framework, particles are approximated by point-particles and the 
forces acting on them are : buoyancy, acceleration of the undisturbed flow, the Stokes drag, 
the added mass and the Basset-Boussinesq history term In this limit an approximate 
equation of motion can be written and coupled to turbulence models or numerical simu- 
lations to predict the motion of the advected inertial particles . This approach 
for instance predicts the wellknown preferential concentration and enhancement of settling 
velocity of heavy particles by the turbulence. Such recent numerical studies Q| have also 
shown to correctly describe acceleration statistics measured for small water droplets (with 
sub-Kolmogorov size) transported in a turbulent air flow ij]. However when the particle 
size is of the order of Kolmogorov scale or larger, such point-particle approximations might 
become inapropriate. 

In this letter, we present an experimental study of the role played by particles finite 
size on their Lagrangian small scale dynamics as they are advected by the turbulence. In 
order to decouple finite size effects from density and seeding concentration effects, we only 
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consider here neutrally buoyant particles (the case of heavy particles will be addressed in a 
forthcoming work) with an extremely low seeding density. Particles are rigid (we use soap 
bubbles, with Weber number much smaller than one in experimental conditions), spherical 
and characterized only by the ratio <fi of their diameter D to the Kolmogorov scale of the 
carrier flow, <p = D/rj. We vary the particles diameter over a range corresponding to inertial 
scales (0 > 1). 

An intuitive phenomenology for finite size effects on the small scales motion of finite 
sized particles advected by a turbulent flow relies on the idea that such particles are not 
sensitive to flow gradients at scales smaller than their own size. In this phenomenology, 
the turbulent energy spectrum of spatial scales actually forcing the particles dynamics is 
therefore somehow filtered, as only turbulent structures at scales typically larger than D do 
contribute to the advection of the particles. We address here the question of the relevance 
of this phenomenology and of how the Langrangian dynamics of the particles (caracterized 
mostly by their acceleration) is affected by such spatial filtering of the turbulent field due 
to the particles finite size. 

Our experiment runs in a large closed wind tunnel with a measurement section of 0.75 mx 
0.75 m where the turbulence is generated behind a grid with a mesh size of 6 cm and 
reproduces almost ideal isotropic turbulence. The results reported here were obtained with 
a mean velocity of the fluid U — 15 m • s _1 and a turbulence level u rms /U ~ 3%. The 
corresponding Reynolds number, based on Taylor microscale, is of the order of R\ = 160. 
The dissipation scale 77 is 240 /im and the energy injection scale L is 6 cm. As particles we 
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FIG. 1: Principle of acoustical Doppler velocimetry : an ultrasonic plane wave of frequency vq is 
emitted in the direction uq by a transducer. A receiver records the wave scattered in the direction 
n by a moving particle with velocity u. The component Uu of the particle velocity along the 
scattering vector q sca t is directly measured by the Doppler shift of the frequency of the recorded 
scattered wave. 



3 



use soap bubbles inflated with helium in order to adjust the density to neutrally buoyant in 
air 6[. The size of the bubbles can be adjusted between 1.5 mm and 6 mm, corresponding 
to ratios to Kolmogorov scale ranging from 7 to 25. For a given setting of the bubble 
size, the standard deviation of the size distribution is relatively low (typically of the order 
of D rms ~ 150 /im), except for the smallest bubbles (D ~ 1.5 mm) for which the nozzle 
generating bubbles has to be pushed to an extreme regime where it becomes slightly unstable 
and produces a wider distribution of particles diameter (typically D 

rms ~ 500 /im in this 

regime). All the different particles we have tested have a diameter which corresponds to 
the carrier turbulence inertial range scales. The smallest particles (0 ~ 7) approach the 
dissipation scale while the largest ones (0 ~ 25) approach the injection scale. The seeding 
density is extremely low (particles are injected individually), so that no global modulation 
of the turbulence by the particles is expected. Particles are individually tracked by ID 
Lagrangian acoustic Doppler velocimetry 0, 0| (figure [[]). We measure the longitudinal 
(streamwise) velocity component v z of the particles along their trajectory. As they travel 
accross the measurement volume, the particles can be tracked during approximately 50 ms 
at most, which corresponds to several dissipation time scales t v = (zz/e) 1 / 2 ~ 3.8 ms [y = 
1.5 • 10 -5 m 2 • s -1 is the air viscosity at working temperature and e ~ 1.0 m 2 • s~ 3 the average 
energy dissipation rate per unit mass at the measurement volume location). For each set of 
particle sizes, we record approximately 4000 tracks at a samplimg rate of 32768 Hz, giving 
more than 10 6 data points for each set. 

In order to characterize the particles dynamics at different time-scales, we analyze the 
Lagrangian velocity increments 5 T v z (t) = v z (t+r)—v z (t). As an example, figured represents 
the probability density functions (PDF) of the Lagrangian velocity increments normalized 
by their standard deviation for increasing values of r for the largest particles (<fr = 25, 
D = 6 mm). The continuous deformation of the PDF from gaussian at large time-scales to 
the development of stretched exponential tails at dissipative time-scales, is the signature of 
an intermittent Lagrangian dynamics (si, [l(| , which we found here to remain present even for 
particles approaching in size the energy injection scale. The wide stretched exponential tails 
of the PDF for the smallest time lags r show that large particles still experience very high 
acceleration events with a probability much higher than gaussian, as previously reported 



for small tracers 



ll|. Here after, we focus our analysis on the acceleration statistics of the 



particles (which is of particular interest since it directly reflects the actual forcing particles 
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FIG. 2: (a) Lagrangian velocity increments probability density functions (PDF) for 6 mm bub- 
bles (<fi ~ 25.1) for different time lags r = r^/50, 2t v , At^, 6t v , Sr^andlOr^ from top to bottom 
; curves have been shifted vertically for readability ; the dashed line represents a gaussian dis- 
tribution of variance 1. (b) Non normalized acceleration PDF for 3 different sizes of particles 
((f) = 7.1; 16.6 and 25.1). (c) Collapse of the acceleration normalized (to variance 1) PDFs for 
different particles sizes (4> = 7.1 — > 25.1) ; the dashed line represents a fit by relation ([1]) with 
s = 0.62. 



experience from the turbulent carrier flow) and on how the turbulent forcing acting on a 
particle depends on its size. Technically, acceleration is obtained by convolution of the 
Doppler velocity signal with a differentiated gaussian kernel. This procedure simultaneously 
differentiates and niters the signal from unavoidable experimental noise 121 ] . 

Figure [2]d represents the acceleration PDF (here acceleration has not been normalized by 
its standard deviation) for three different particle sizes. Acceleration PDF clearly narrows 
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and peaks around zero as the particle dimension increases. This gives a clear evidence that 
depending on their size, particles experience a turbulent forcing with different statistical sig- 
nature. Larger particles present a narrower acceleration distribution than smaller, meaning 
that they are less sensitive to extreme turbulent events. An interesting observation is that 
the evolution of the PDF with the particle diameter is mostly associated solely to a decrease 
of the acceleration variance of the particles as their diameter increases. Indeed, when we 
normalize the acceleration component a z of the particles by its standard deviation a z>rms , the 
corresponding PDFs are found to collapse for all particle diameters (figure [2b). Therefore, 
the global shape of the acceleration PDF, normalized to variance 1, is found to be preserved 
for all particle sizes. (We note that the PDF for the smallest particles (0 = 7.1) does not 
collapse as well as the others, especially the PDF tails deviate. This is very likely due to the 
larger dispersion in the size distribution of the smallest particles which tends to overstimate 
the actual acceleration standard deviation and to alter the global shape of the PDF.) Over 
the range \a z \ < 10 a 2jrms the shape of the normalized PDF is well described by the function 



p 3s 2 /2 

V(x) = — =- 
4\/3 
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V2s 



(1) 



where s is a parameter related to the global shape of the PDF. This form for the PDF of 
an acceleration component was previously suggested for the case of fluid tracer particles 
(with size comparable to r/) [3| and was shown to be related to a lognormal distribution 
(of variance s 2 ) of the whole vector acceleration magnitude [ijj] under the condition of 
isotropic turbulence. We found here that this form still correctly describes the acceleration 
component statistics of finite sized particles and considering the good isotropic conditions 
in our experiment, we can expect acceleration magnitude of finite sized particles to be 
indeed well described by a lognormal distribution. Our experimental data is best fitted 
with a value for the parameter s in relation ([T|) equal to 0.62, which corresponds to a 
flatness for the acceleration fluctuations equal to 8.4 (in good agreement with the direct 
estimation of the flatness from the data). This value is significantly smaller than values 

nn 

previously reported for fluid tracers (of the order of 20) at similar Reynolds number |9|, 115] . 
although a unambiguous comparison between those experiments and ours is made difficult 
due experimental specificities ; in particular while our grid generated turbulence is largely 
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isotropic, the von Karman flow studied in [9|, [15| has a large scale anisotropy, which have been 



shown to possibly affect small turbulent scales, even in high Reynolds numbers regimes [16]. 
If, in spite of these differences, a continuity between those previous experiments and ours 
is however attempted, the difference between the acceleration flatness between small fluid 
tracers and our finitie sized particles would mean that a rapid decrease of the acceleration 
flatness must occur for particles sizes around a few 77 (0 ~ (9(1)) as the size increases from 
dissipative to inertial range scales. In any case, our results show that acceleration statistics 
reach an asymptotic regime for inertial range sizes, regardless of the particles size (at least 
up to ~ 30), described by the PDF in figure [2b and equation (CD). An interesting and new 
result to be stressed is that the same identical non-gaussian acceleration PDF is preserved 
even for particles with size approaching injection scales. Based on the classical image of 
Eulerian intermittency of the carrier turbulent velocity field (where velocity increments are 
non-gaussian at small scales and tend continuously toward a gaussian at larger scales), one 
might have accordingly expected a monotical trend of particles acceleration statistics toward 
a gaussian as particles size approaches large scales. Similarly, one might have expected the 
acceleration statistics of particles of size D to be somehow related to the Lagrangian velocity 
increments statistics calculated along fluid tracers trajectories at a time lag given by the 



eddy turnover time at scale D, which because of Lagrangian velocity intermittency [9|, Il0j . 
would also tend continuously toward a gaussian at large scales. This not being the case, 
we can conclude that the turbulent forcing of finite sized particles with inertial range sizes 
is not trivially related to intermittency (neither Eulerian nor Lagrangian) of the turbulent 
velocity field around the particle; in particular acceleration PDF for particles with a given 
size cannot be simply deduced from statistics of velocity increments at scales corresponding 
to the particle size. 

A possible reason to this is that from a dynamical point of vue, acceleration reflects 
the forcing felt by the particle, which, for inertial range scales, is mostly due to turbulent 
fluctuations of the pressure gradients in the vicinity of the particule. Under an ergodic 
assumption, particles acceleration statistics can therefore be related to Eulerian pressure 
increments statistics, which are indeed not trivially related to velocity increments statistics 
\v\ . In particular, in this context the acceleration variance of a particle with diameter D is 
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related to the second moment of the pressure increments at the scale of the particle: 

<«2WJ»«^. ( 2 ) 

where 5"^(r) = ((p(x + r ) — p(x)) 2 ), is the Eulerian second order pressure structure function 
(which only depends on r = \r\ in homogeneous isotropic conditions). In the frame of 
Kolmogorov 41 phenomenology, the classical inertial scaling for Sf (r) is (er) 4 / 3 
(or equivalently e 4//3 A;~ 7 / 3 for the pressure spectrum in k wave number space). Although 
this scaling is still controversial Q, 0], there are experimental 21, 22] and numerical [20I ] 
evidences suggesting that it is most likely correct in ranges of Reynolds numbers and scales 
consistent with our flow and particle sizes. As a consequence, the acceleration variance of 
particles with diameter D in inertial range scales should follow the scaling: 

k>^m.(o) = «; £ 1/3 b- 2/s . p) 



particle 

This scaling can also be directly derived by simp 



e dimensional considerations in the frame- 



work of Kolmogorov's inertial range hypotheses 1^] , but we stress here the dynamical cor- 
respondance with the pressure increments statistics. Relation (j3J) predicts a decrease of the 
acceleration variance with increasing particles diameter. Our measurements show indeed 
that while the global shape of the normalized acceleration PDF is found not to depend on 
particle size, a 40% decrease is measured for the acceleration variance itself between the 
smallest and the largest particles we have tested. Figure [3h shows more precisely the de- 
crease of the acceleration variance (made dimensionless by e 4 / 3 // -2 / 3 ) with the particle size. 
Errorbars are mostly associated to the particles size distribution and to uncertainties in the 
determination of e. The decay of acceleration variance is found in good agreement with 
relation (jHJ) for particles larger than about 4 mm (<fr > 15), as it can be more clearly seen 
on the compensated plot on figure [3b, where an asymptotic plateau (around a' ~ 18) is 
reached for the largest particles. 

In the limit of very small particles (ideally <p <C 1) acceleration variance should saturate 
to an asymptotic value corresponding to the intrinsinc acceleration variance of the turbulent 
flow itself. In this limit, scaling ([3]) which has been derived from inertial range considera- 
tions, should be replaced by its equivalent small scale (dissipative) form, given by the usual 
Heiseinberg-Yaglom relation for fluid particles [ijj]: 
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$ = D/ri t|> = D / rj 

FIG. 3: (a) Dimensionless acceleration variance as a function of particle size ; errorbars are mostly 
associated to the particle size distribution and to the uncertainty on the determination of e. The 
point for the smallest value of has larger errorbars due bubble nozzle instability which produces 
a larger size distribution in this regime, tending very likely to overestimate particles acceleration 
variance. (b) Compensated acceleration variance as a function of particle size. In both figures the 
solid line represents a D~ 2 ^ 3 law. 



<«•>„ = a„ £ "/V^. (4) 

Using high speed optical particle tracking Voth and collaborators [l5| have shown that accel- 
eration variance of particles such that < 5 is already correctly described by the dissipative 
scaling (j3J) with a Reynolds number dependent constant a® (reasonably, the dissipative scal- 
ing is reached far before the ideal limit <C 1). Figure shows a clear deviation from 
the inertial scaling (J3J) as soon as < 15- Therefore, we can conclude that the crossover 
between the inertial scaling (j3J) and the dissipative scaling (j3j) for the acceleration variance 
occurs for particles with size in the range 5 < < 15. Although we are not able here to 
study particles smaller than < 7 (for technical reasons related to our bubble generator 
nozzle), an extrapolation to the limit 0-^0 from figure [3^ gives a rough estimate for the 
constant a in (J2J) at the working Reynolds number (R\ ~ 160) of the order of On~ 3 which 
is in good agreement with the experimental value by Voth and collaborators 15j and from 
earlier predictions from direct numerical simulations [19I. I20I] . 

We have studied Lagrangian statistics of finite sized neutrally buoyant particles trans- 
ported in a trubulent flow of air, with a particular focus on particles acceleration. These 
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are the first systematic experimental measurements of finite size effects on the transport of 
material particles ever done in an open and nearly isotropic turbulent flow ; they were made 
possible by the combination of a novel acoustical Lagrangian tracking system and a versatile 
particle generator. We have shown that particle dynamics displays intermittency and that 
particle acceleration statistics is described by the same probability density function for all 
particles with inertial range sizes, where only the acceleration variance displays a significant 
dependence on particle size. In particular acceleration fluctuations remain non-gaussian 
even for large particles (approaching the injection scale). The decrease of the acceleration 
variance with particle size is in good agreement with classical inertial scaling for pressure 
increments in the carrier flow, what supports the connection between the Lagrangian dy- 
namics of the particles and Eulerian statistics of the carrier turbulent flow via the pressure 
increments field. A deeper theoretical and numerical insight on pressure increments statis- 
tics at different spatial scales (for instance in terms of a systematic analysis of pressure 
increments PDF) would very likely help the further understanding and the modeling of such 
finite size effects. Our results give an experimental support to the idea, generally assumed 
for the modelling of particles transport, that turbulent pressure fluctuations only at scales 
larger than the particles size do contribute to the transport and the dispersion of the par- 
ticles, which is of particular interest in the frame of Large Eddy Simulations (LES) for the 
study of particles-turbulence interactions. Finally, in a purely hydrodynamics context, the 
present experiments also open new possibilities for characterizing Eulerian statistics of tur- 
bulent pressure fluctuations at inertial scales by means of Lagrangian measurements of the 
acceleration of finite sized particles. 
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